Test the time-reversal modified universality of the Sivers function 
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We derive the time-reversal modified universality for both quark and gluon Sivers function from 
the parity and time-reversal invariance of QCD. We calculate the single transverse-spin asymmetry 
of inclusive lepton from the decay of W bosons in polarized proton-proton collision at RHIC, in terms 
of the Sivers function. We find that although the asymmetry is diluted from the W decay, the lepton 
asymmetry is at the level of several percent and is measurable for a good range of lepton rapidity 
at RHIC. We argue that this measurable lepton asymmetry at RHIC is an excellent observable for 
testing the time-reversal modified universality of the Sivers function. 
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I. Introduction. Much of the predictive power of 
perturbative Quantum Chromodynamics (QCD) is con- 
tained in factorization theorems They normally in- 
clude two assertions. One is that a physically measured 
quantity can be factorized into some perturbatively cal- 
culable short-distance hard parts convoluted with non- 
perturbative long-distance distribution functions (or ma- 
trix elements) . The other is the universality of the non- 
perturbative distribution functions. Predictions follow 
when processes with different hard scatterings but the 
same distribution functions are compared. With one set 
of universal parton distribution functions (PDFs) the 
leading power collinear QCD factorization formalisms 
have been very successful in interpreting and predicting 
almost all existing data from high energy collisions with 
momentum transfers larger than a few GcV [2] . 

The phenomenon of single transverse-spin asymmetry 
(SSA), A N = (<t(Sl) - o-{-S ± ))/(a(S x ) + a(-S ± )), de- 
fined as the ratio of the difference and the sum of the 
cross sections when the single transverse spin vector S± is 
flipped, was first observed in the hadronic A production 
at Fermilab in 1976 as a surprise [3J. Large SSAs, as large 
as 30 percent, have been consistently observed in various 
experiments involving one polarized hadron at different 
collision energies jj]. The size of the SSAs presented a 
challenge to the leading power collinear QCD factoriza- 
tion formalism 5] , and provided an excellent opportunity 
to probe a new domain of QCD dynamics. 

Two widely discussed theoretical approaches have been 
proposed to evaluate the observed SSAs in QCD. One 
generalizes the QCD collinear factorization approach to 
the next-to- leading power in the momentum transfer Q , 
and attributes the SSA to the quantum interference of 
scattering amplitudes with different numbers of active 
partons [3, [!]• The size of the asymmetry is deter- 
mined by new three-parton correlation functions . This 
generalized collinear factorization approach is more rel- 
evant for the SSAs of cross sections whose momentum 
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transfers Q ^> Aqcd- The other approach factorizes 
the cross sections so as the SSAs in terms of the trans- 
verse momentum dependent (TMD) parton distributions 
[10, [H El El El, and attributes the SSAs to the 
Sivers function [l5l | (or the Collins function if a final- 
state hadron was observed 16]). The TMD factoriza- 
tion approach is more suitable for the SSAs of cross sec- 
tions with two very different momentum transfer scales, 
Qi ^S> Q2 Aqcd- These two approaches each have 
their kinematic domain of validity, they were shown to 
be consistent with each other in the kinematic regime 
where they both apply (l7| . 

However, there is one crucial difference between these 
two factorization approaches besides the difference in 
kinematic regimes where they apply. The Sivers function 
in the TMD factorization approach could be process de- 
pendent, while all distribution functions in the collinear 
factorization approach are universal. It was predicted 
by Collins [10] on the basis of time-reversal arguments 
that the quark Sivers function in semi-inclusive deep in- 
elastic scattering (SIDIS) and in Drell-Yan process (DY) 
have the same functional form but an opposite sign, a 
time-reversal modified universality. In this Letter, we 
derive the same time-reversal modified universality for 
both quark and gluon Sivers function from the parity 
and time-reversal invariance of QCD. 

The experimental check of this time-reversal modified 
universality of the Sivers function would provide a criti- 
cal test for the TMD factorization approach [H EH, El 
E3I EH- Recently, the quark Sivers function has been 
extracted from data of SIDIS experiments [Hj]. Future 
measurements of the SSAs in DY production have been 
planned [h|. In this Letter, we present our calculation of 
the SSAs of inclusive single lepton production from the 
decay of W bosons. The W production and DY share the 
same Sivers function. We find that although the asym- 
metry is diluted from the decay of W bosons, the lepton 
asymmetry is significant and measurable for a good range 
of lepton rapidity at RHIC. We show that the lepton 
SSAs provide the better flavor separation of the quark 
Sivers function than what the standard DY can do. We 
also show that the lepton SSAs are sharply peaked at 
transverse momentum px ~ Mw/2 with W mass Myy. 
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Since leptons from heavy quarkonium decay and other 
potential backgrounds are unlikely to be peaked at the 
Pt ~ Mw/2, we argue that the SSA of inclusive high pt 
leptons at RHIC is an excellent observable for testing the 
time-reversal modified universality of the Sivers function. 

II. The QCD prediction. The predictive power of the 
TMD factorization approach to the SSAs relies on the 
universality of the TMD parton distributions. For the 
lepton-hadron SIDIS, t(l) + h(p, S) -> i'(V) + h'{p') + X, 
the factorized TMD quark distribution has the following 
gauge invariant operator definition [2fJ, 



f^l s (x,k ± ,S) = 



dy d 2 y ± Jxp + y -- i]i± . yj 



(2tt) 3 

x (p,^(0-,0 ± )*t({oo,0},Oj 



< (oo,{yx,Ox» 



+ 



(1) 



x $ n ({oo,y },yx)ip(y ,y±)\p,S), 

where y + = + dependence is suppressed and the gauge 
links from the final-state interaction of SIDIS are 



$„({oo,y },y x ) = Ve 
$ ni (oc,{yx,Ox}) = Ve 



-ig dy 1 A fl (y 1 ,y ± ) 

, 



where V indicates the path ordering and the direction 
nx is pointed from 0j_ to yj_. Here we define the light- 
cone vectors, = (n + ,n~,nx) = (0,1, Ox) and = 
(1,0, Ox), which project out the light-cone components 
of any four- vector V^ 1 as V ■ n — V + and V ■ n = V~. 

For the DY, h(p,S) + h'{p') -> 7*(Q)[-» + X, 

the factorized TMD quark distribution is given by 



/.j/Xt ( x , kx, S) 



dy d 2 y± 



-ik_L-yj 



(27t) 3 

x (p ) 5|^(0-,0x)4({-oo,0},0 J 
x $t_ L (_ cx >,{y±,0 ± })^- 



(3) 



x $„({-oo,y },yx)VXy ,y±)\p,S) 

where the past pointing gauge links were caused by the 
initial-state interactions of DY production From 
Eqs. H]) and ©, it is easy to show that the collinear 
quark distributions are process independent, 



d 2 k ± f^l s (x,k ± ,S) = d 2 k ± f^ hr (x,k ± ,S), (4) 



if the same renormalization scheme was used for the ul- 
traviolet divergence of the kx integration. 

Let | a) = \p, S) and (/3| be equal to the rest of the 
matrix element in Eq. ([T]) 9]. From the parity and 
time-reversal invariance of QCD, (ap\(3p) — (a|/3) and 
(Pt\olt) — ((x\P), where \ap) and |/3p), and \o>t} and 
\Pt) are the parity and time-reversal transformed states 
from the states \a) and respectively, we derive 



f^l s (x, kx, S) = /°* t (x, kx, -S) . 



and conclude that the spin-averaged TMD quark distri- 
butions are process independent. Following the notation 
of Ref. [3 , we expand the TMD quark distribution as 



f q /ht{x,kx,S) = f q / h (x,k ± ) 



±A N f q/h r(x,k ± )§-(pxt ± ) (6) 



where k± = |kx|, p and kx are the unit vectors of p and 
kx, respectively, f q /h(x, k±) is the spin-averaged TMD 
distribution, and A N f q / h r (x, k±) is the Sivers function 
[H[. Substituting Eq. J6) into Eq. ©, we obtain, 



(7) 



which confirms the Collins' prediction [10] that the Sivers 
function in SIDIS and in DY differ by a sign. 

We define the gauge invariant TMD gluon distribu- 
tion in SIDIS and in DY by replacing the quark oper- 
ator in Eqs. fT} and §5§ by the gluon opera- 
tor F +ll F +1/ {—g tlv ), and the gauge links by those in the 
adjoint representation of SU(3) color. From the parity 
and time-reversal invariance of the matrix elements of 
the TMD gluon distribution, we find, like Eq. ([5]), 



fg}h1 S ( X ' kx, S) 



fg/h.1 i X ' kx, —S) 



(8) 



Applying Eq. © to the gluon TMD distribution, we de- 
rive the same time-reversal modified universality for the 
gluon Sivers function, 



(9) 



The sign change of the Sivers function is a property of 
the gauge invariant TMD parton distributions. 

777. Lepton SSAs from W production. The SSAs of W 
production at RHIC, A(p A ,S±) + B(p B ) -> W ± {q) -> 
^ ± (p) + A, were proposed in Refs. [2l| to measure the 
Sivers function. However, it is difficult to reconstruct 
W bosons by the current detectors at RHIC. We present 
here our predictions for the SSAs of inclusive lepton pro- 
duction from the decay of W bosons. We also present 
the SSAs of W production for a comparison. 

We use the TMD factorization formalism because W 
bosons at RHIC are likely produced with transverse mo- 
mentum |qx| -C M\y- We work in a frame in which the 
polarized hadron A moves in the -l-z-direction. We have 
the leading order spin-averaged W cross section 



da A 



B-*W 



dy w d 2 q± 



<r J2\ Vab \ 2 / d 2 k a ^d 2 k b± f a/A { 



xfb/B{xb,h±_)5 2 (ci A 



), (10) 



where yw is the W rapidity, o- — (n/3) \/2Gf Myy/s is 
the lowest order partonic cross section with the Fermi 
weak coupling constant and s = (p A +Ps) 2 , J2ab 
(5) runs over all light (anti)quark flavors, V a b are the CKM 
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matrix elements for the weak interaction. The parton 
momentum fractions in Eq. (jlOp are given by 



M 



M 



Xb 



w e -vw 



(11) 



to the leading power in q^/M^. Similarly, we have the 
leading order factorized spin-dependent W cross section 
Aa(S±) = [a(S±) - a(-S±)]/2 as 



dAtr A TB->w(S±) _ &o_ 
dy w d 2 q± 2 



^2\v ab \ 2 fd 2 k a± d 2 k b 



x S ± ■ (p A x k a± ) A N fVf At (x a , k a± ) 
x fb/B(xb, k b ±)5 2 (<i± - k a ± - k b ±) . (12) 
The SSA of W production is then defined as, 

dAa(S±_)AT b^w / da AB ^ w 



A 



(W) 
N 



dy w d 2 qj 



dy w cPqj 



(13) 



whose sign depends on the sign of the Sivers function and 
the direction of the spin vector S±. 

To evaluate the SSA in Eq. (fT3")l , we use the parame- 
terization of TMD parton distributions in Ref. [18| , 



f q /h(x,k±) = f q (x) 



1 



-kl/(k 2 ± ) 



(14) 



h(k±_) 



n(kl) 

2JV q (x)h(k x )f q/h (x,k x ), (15) 



Mi 



(16) 



A™ = S ± ■ (p A x qj.) 



i(k 2 



<*1 >-<*=;> 



c=1 >+<*s> 



2e Ea 6 1^6 



T- (fc s 2 )] 2 

Na(Xa)} fa(x a ) fb(Xb) 



Mi 



Eafcl^b| 2 /a(^a)/b(a; & ) 



2<fc1 > 



(17) 



where (A; 2 ) = Mj 2 (fcj_)/[Mf + (fc^)] and the "-" sign in 
front of Af a (x a ) is from Eq. ([7]). If we choose the S± along 
the y-axis as in Ref. 5j_-(p^xqj_) = qTCOs(4>w) with 
<7t = |q_i_| and azimuthal angle 4>w- For our numerical 
predictions below, we choose <f>w = and the GRV98LO 
parton distribution [22] for f q (x) to be consistent with 
the usage of the TMD distributions of Ref. [H]. 

In Figs. [1] and [21 we plot the A N from Eq. (JT7J) 
at </s — 500 GeV. The W asymmetry is peaked at 
or <C Mw and is much larger than that of DY production 
[19(. This is because the u and d Sivers functions have 
an opposite sign, and they partially cancel each other in 
their contribution to the DY asymmetry, while they con- 
tribute to the W + and W~ separately. The large W~ 
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FIG. 1: An as a function of IF-boson rapidity. 
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where f q (x) is the standard unpolarized parton distribu- 
tion of flavor q, (k 2 ) and M\ are fitting parameters, and 
Af q (x) is a fitted distribution given in Ref. [HI]. By car- 
rying out the integration G? 2 k a ^d 2 k&^ in Eqs. (flT)]) and da 
(fi"2"j) analytically, we obtain, 



FIG. 2: An as a function of IF-boson transverse momentum. 

asymmetry is caused by a large d Sivers function [l8j |. 
The negative d Sivers function in SIDIS gives the posi- 
tive W~ asymmetry. The rapidity dependence in Fig. [T] 
provides excellent informations for the flavor separation 
as well as the functional form of the Sivers function if we 
could reconstruct the W bosons. 

After integrating over the momentum of (anti)ncutrino 
from the W decay, we obtain the leading order factorized 
cross section for the production of leptons of rapidity y 
and transverse momentum pj_, 



ATB^£(p) 



(S±) 



dyd 2 pA 



E 

a . b 
pDY 



V ab 



dx a d 2 k al _ / dx b d 2 k b ± 



X / a /AT( a; a,k ai ,5 i ) fb/B(xb,h±) 
1 . |2 

Mab^A S(s + t + u), (18) 



167T 2 J 



where s, t, and u are the Mandelstam vaeriables and 
the leading order partonic scattering amplitude square, 
|A4 Q i,_.e| , is given by 



8(G F M W ) 2 



(s-M w r+M w T 2 w 



(19) 



for partonic channels ab = du, su, du, su; or by the same 
one with the u 2 replaced by P for the rest light flavor 
channels ab = ud, us, ud, us. Tw in Eq. (|19|) is the W lep- 
tonic decay width. Substituting Eq. ((6]) into Eq. (fl8|) . we 
derive both the spin-averaged and spin-dependent cross 
sections, from which we evaluate the SSAs of inclusive 
lepton production from W decay numerically. 

In Figs. [3] and 31 we present our predictions for the in- 
clusive lepton asymmetry from the decay of W bosons at 
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FIG. 3: An as a function of lepton rapidity. 
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FIG. 4: An as a function of lepton transverse momentum. 

RHIC energy. Although the decay diluted the size of the 
asymmetry, the lepton inherited all key features of the W 
asymmetry in Figs. [1] and O As shown in Fig.rjJ the lep- 
ton asymmetry is sharply peaked at px ~ 41 GeV, which 
should help control the potential background. The differ- 
ence in rapidity dependence of the W + and W~ in Fig. [3] 
provides the excellent flavor separation of the Sivers func- 



tion, as well as rich information on the functional form. 
For a good range of rapidity, the lepton asymmetry is 
measurable at RHIC. 

IV. Summary and Conclusions. In summary, we have 
derived the time-reversal modified universality for both 
quark and gluon Sivers functions from the parity and 
time-reversal invariance of the gauge invariant matrix el- 
ements that define the TMD parton distributions. We 
confirm the Collins' prediction for the sign change of the 
quark Sivers function in SIDIS and in DY [n|. The sign 
change of the Sivers function in SIDIS and in DY is a nat- 
ural property of the gauge invariant TMD parton distri- 
butions in QCD. Corresponding sign change of the SSAs, 
if they could be factorized in terms of these TMD parton 
distributions, is a fundamental prediction of QCD. 

We have calculated, in terms of the TMD parton dis- 
tributions, the SSAs of W production as well as inclu- 
sive lepton production from the decay of W bosons in 
polarized proton-proton collision at RHIC energy. We 
find that although the asymmetry is diluted from the W 
decay, the lepton asymmetry is at the level of several 
percent and measurable for a good range of lepton rapid- 
ity at RHIC. Because the lepton asymmetry is sharply 
peaked at the px ~ 41 GeV, the potential background 
could be strongly suppressed. We conclude that this 
measurable lepton asymmetry at high px at RHIC is an 
excellent observable for measuring the Sivers functions of 
different flavors and for testing the time-reversal modified 
universality of the Sivers function. 
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[1] for reviews, see: J. C. Collins, D. E. Soper and G. Ster- 

man, Adv. Ser. Direct. High Energy Phys. 5, 1 (1988). 
[2] for reviews, see: J. Pumplin, D. R. Stump, J. Huston, 

H. L. Lai, P. Nadolsky and W. K. Tung, JHEP 0207, 

012 (2002); A. D. Martin, R. G. Roberts, W. J. Stirling 

and R. S. Thorne, Eur. Phys. J. C 23, 73 (2002). 
[3] G. Bunce et ai, Phys. Rev. Lett. 36, 1113 (1976). 
[4] for reviews, see: U. D'Alesio and F. Murgia, Prog. Part. 

Nucl. Phys. 61, 394 (2008). 
[5] G. L. Kane, J. Pumplin and W. Repko, Phys. Rev. Lett. 

41, 1689 (1978). 
[6] J. W. Qiu and G. Sterman, AIP Conf. Proc. 223, 249 

(1991); Nucl. Phys. B 353, 137 (1991). 
[7] A. V. Efremov and O. V. Teryaev, Sov. J. Nucl. Phys. 

36, 140 (1982). Phys. Lett. B 150, 383 (1985). 
[8] J. W. Qiu and G. Sterman, Phys. Rev. Lett. 67, 2264 

(1991), Nucl. Phys. B 378, 52 (1992). 
[9] for example, see: Z. B. Kang and J. W. Qiu, Phys. Rev. 

D 79, 016003 (2009), and references therein. 
[10] J. C. Collins, Phys. Lett. B 536, 43 (2002). 
[11] D. Boer, P. J. Mulders and F. Pijlman, Nucl. Phys. B 

667, 201 (2003). 



[12] X. d. Ji, J. P. Ma and F. Yuan, Phys. Lett. B 597, 299 

(2004), Phys. Rev. D 71, 034005 (2005). 
[13] J. C. Collins and A. Metz, Phys. Rev. Lett. 93, 252001 

(2004). 

[14] J. Collins and J. W. Qiu, Phys. Rev. D 75, 114014 (2007). 

[15] D. W. Sivers, Phys. Rev. D 41, 83 (1990), 43, 261 (1991). 

[16] J. C. Collins, Nucl. Phys. B 396, 161 (1993). 

[17] X. Ji, J. W. Qiu, W. Vogelsang and F. Yuan, Phys. Rev. 
Lett. 97, 082002 (2006), Phys. Rev. D 73, 094017 (2006), 
Phys. Lett. B 638, 178 (2006); Y. Koike, W. Vogelsang 
and F. Yuan, Phys. Lett. B 659, 878 (2008). 

[18] M. Anselmino et ai, Eur. Phys. J. A 39, 89 (2009) 

[19] J. C. Collins et ai, Phys. Rev. D 73, 094023 (2006). 

[20] X. d. Ji and F. Yuan, Phys. Lett. B 543, 66 (2002); 
A. V. Belitsky, X. Ji and F. Yuan, Nucl. Phys. B 656, 
165 (2003). 

[21] S. J. Brodsky, D. S. Hwang and I. Schmidt, Phys. Lett. 

B 553, 223 (2003); I. Schmidt and J. Soffer, Phys. Lett. 

B 563, 179 (2003). 
[22] M. Gluck, E. Reya and A. Vogt, Eur. Phys. J. C 5, 461 

(1998). 



